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Abstract 

A statistical analysis of the spectrum of two particle - two hole doorway states in 
a finite nucleus is performed. On the unperturbed mean-field level sizable attractive 
correlations are present in such a spectrum. Including particle-hole rescattering ef- 
fects via the residual interaction introduces repulsive dynamical correlations which 
generate the fluctuation properties characteristic of the Gaussian Orthogonal En- 
semble. This signals that the underlying dynamics becomes chaotic. This feature 
turns out to be independent of the detailed form of the residual interaction and 
hence reflects the generic nature of the fluctuations studied. 
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Theoretical models aiming at the quantal description of an excitation and sub- 
sequent decay of a collective state in a many-body system are usually based on a 
division of the full Hilbert space into two sectors, Si and S 2 , spanned by the vec- 
tors |1) and 1 2), respectively. At the same time, the hamiltonian is represented as 
H = H + V, such that (1\H \1') = e°hi>, <2|J3" |2') = e%6 22 > and (l|J2o|2) = 0. 
When V is taken into account the above relations no longer hold. By diagonalizing 
H in the basis |2) and redefining |2) one can still have (2\H\2') = €2622', however. 

A collective state |/} such as a plasmon or a nuclear vibrational mode is defined 
as an eigenstate of H in the subspace Si: |/) = I n the full space, in- 

cluding Si and S2, I/) is, however, no longer an eigenstate of H but rather a wave 
packet which begins to 'leak' into the space S 2 . This constitutes a mechanism for 
dissipation. Thus, for the time-dependent state one has 

!/(*)> = E /i(*)|l> + E /a(*)|2>, Mt = 0) = 0) (1) 

1 2 

and the Schrodinger equation for f\ and f 2 reads: 

d 
1— 
dt 

Applying a procedure similar to the Nakajima-Zwanzig projection technique O, 
i.e., solving the second of those equations for f 2 (t), inserting into the first one for 
/i(i) and assuming that the basis |2) is already defined such that {2\H\2') = €2622', 
yields 

^/l(*)-E i? H'/l'(*)=E f dTf v (t-r)v lv {r), (3) 
dt y Jo 

where 

v\v (r) = -i E H i2H 2 v exp(-ie 2 r) . (4) 
2 

The r.h.s. in eq. (||) is a generalized non-local collision term which compensates for 
the formal elimination of the space S2. The matrix elements H\ 2 = (1\H\2) describe 
the degree of mixing between Si and S2 after inclusion of V. The internal structure 
of space S2 is represented by the eigenvalues €2 of the full hamiltonian H in this 
space. 

In nuclear physics the space Si is typically spanned by one particle - one hole 
(lplh) states generated from the mean field hamiltonian, H a . Diagonalizing the 
full nuclear hamiltonian in such a basis gives the RPA boson excitations. Ideally, 
the space S2, absorbing such an excitation, should contain all possible npnh states. 
Then, after diagonalization of H, the space S2 would reflect the entire complexity of 
the spectrum of a compound nucleus. This, however, is neither possible for practical 
reasons (too many states) nor does it seem necessary for physical reasons. The 
nuclear hamiltonian involves predominantly two-body interactions and thus couples 
the lplh RPA bosons chiefly to the 2p2h components in the space S 2 - In this way, 
the 2p2h states play a special role of 'doorway states' Q. For these reasons, the 
space S2 is usually specified in terms of the 2p2h states alone Q. The resulting 
methods are known as extended RPA and are believed to account for the spreading 
width of collective modes. Moreover, since the spectrum of such states is already 
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very dense in energy, the interaction in the 2p2h - space is usually neglected which 
corresponds to using the unperturbed 2p2h energies e 2 as e 2 in eq. (||). 

At this point one should recall that 'nuclear dissipation is a process occurring 
in a closed system, and very likely a consequence of chaotic motion' Evidence 
supporting the connection between dissipation and chaos comes both, from classi- 
cal H and from quantum Q considerations. In light of the introductory remarks 
such a link seems indeed to be present. The fluctuation properties of the spectra of 
compound nuclei are consistent with those of the Gaussian Orthogonal Ensemble 
(GOE) of random matrices ||. The same fluctuation properties are identified theo- 
retically H as well as experimentally |T(J for those quantum systems whose classical 
counterparts are chaotic. Similar conclusions can be drawn from the study of open 
phase space, scattering phenomena |]ll[ . It is thus natural to require that modeling 
the space S 2 should preserve this fundamental property of the compound nucleus. 
Because of the special role played by the 2p2h states, the purpose of the present let- 
ter is to quantitatively explore the problem of whether and under which conditions 
already the basis of such states can support the GOE fluctuation characteristics and 
thus makes the underlying dynamics chaotic. 

To make our analysis as realistic as possible the mean field hamiltonian H Q , 
generating the single-particle states, is specified in terms of a local Woods-Saxon 



potential [12]: 



U(r) = -— \^ ~ VL8~(-^^)ff- a) + V c (r) (5) 
1 + e a I u>r 1 + e a 

where V c (r) is the Coulomb potential. The parameters V Q , Vls and o are adjusted 
such as to reproduce the ground state properties of the nucleus and, in particular, 
the single-particle energies near the Fermi surface. To determine the single-particle 
wave functions from this potential, H Q is diagonalized in a harmonic oscillator basis 
including 16 major shells. 

As a residual interaction V we adopt the zero-range Landau-Migdal interaction 

V(ri,r 2 ) = C (/ + f'n ■ t 2 + gai ■ a 2 + 9 <?\ • o" 2 ti • r 2 )<5(n - r 2 ) (6) 

which provides a good description of low-energy nuclear excitations. In extended 
RPA calculations with this interaction (including unperturbed 2p2h states only) 
the basis of single-particle states is typically specified [13] by four major shells, two 



below and two above the Fermi surface. Consistently, in the present study we use 
the basis of similar size. We employ the set of empirical Landau-Migdal parameters 
as given in ref.@ (/ = -0.1,/' = 0.6,5 = 0.15 and g' = 0.7 with C D = 392MeV 
fm 3 ) and perform the calculations for 40 Ca. 

In spherical nuclei, such as 40 Ca, the total angular momentum J and parity 
7T are good quantum numbers which determines the selection of the basic 2p2h 
states, active in the decay of a given RPA boson excitation. All three types of 
the matrix elements generated by the residual two-body interaction in the basis 
of 2p2h states are shown diagrammatically in Fig. 1 and the angular momentum 
coupling scheme is also indicated. The single-particle basis of four major shells 
in the Woods-Saxon potential allows realistic estimates (l3| of global spreading 
widths because the number of the corresponding 2p2h states is already of the order 
of 10 3 - 10 4 , depending on the multipolarity of an excitation with the maximum 
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of 11720 for J n = 3 _ . The matrix of such a size cannot, however, be diagonalized 
with satisfactory precision. A manageable number (2696) of 2p2h states is found for 
J 71 " = + and, therefore, our further discussion will be limited to this multipole. The 
states are not coupled to the good isospin T because the hamiltonian, as specified 
above, explicitly violates the isospin symmetry. An even weak violation of this 
symmetry is known [|lf| to correlate the T = and T = 1 states which allows to 
combine the relevant spectra for a statistical analysis. 

The most obvious measure of spectral fluctuations is the one expressed in terms 
of the nearest-neighbor spacing (NNS) distribution. A standard procedure of anal- 
ysis is to normalize the spectrum such that the fluctuations on different energy 
scales are directly comparable. Here we perform the corresponding unfolding 16] 
by approximating the integrated density of states with polynomials up to the order 
of 12. This guarantees stability of the result. Before unfolding we discard the 50 
lowest states in order to eliminate those states which are located far from their main 
concentration. Similarly we discard the 50 highest-lying states. 

We find three qualitatively different situations as illustrated in Fig. 2. Part (a) 
shows the NNS distribution for the spectrum of unperturbed 2p2h states coupled 
to J n = + . Such a spectrum is not generic and is characteristic for a narrow class 
of integrable systems involving extra correlations in the hamiltonian [17|. The pro- 
nounced peak for small nearest-neighbor separations illustrates a strong tendency 
of states for clustering. Actually, in most cases these are even exact degeneracies. 
They reflect the fact that within the single-particle basis used there are many more 
two hole (hh) states than the number of different energies available. Analogous 
partitioning is even more restrictive on the two particle (pp) side. Including the 
interaction in the pp (diagram (a) in Fig. 1) and hh (diagram (c)) channels removes 
those degeneracies and, as shown in part (b) of Fig. 2, immediately brings the spec- 
trum to the known universality class of generic integrable systems jL7j] characterized 
by a completely uncorrelated sequence of eigenenergies. As a consequence, the NNS 
follows a Poisson distribution. This result essentially does not depend on precise 
values of the parameters of the residual interaction. Already, a comparatively weak 
perturbation (for instance 20% of the original strength) produces a similar picture. 
The other important element which makes the spectrum uncorrelated is that, at this 
level, the system considered (2p2h) still remains a simple product of two subsystems 
(pp and hh) and the corresponding energies are the sums of two independent com- 
ponents. The crucial step which brings us to the situation displayed in part (c) of 
Fig. 2 is the inclusion of the ph-type matrix elements, represented by the diagram 
(b) of Fig. 1. These matrix elements introduce such important correlations that 
the resulting NNS distribution almost perfectly follows a Wigner distribution, i.e. 
is consistent with the GOE. In fact, including only this diagram and ignoring the 
other two ((a) and (c)) gives the same result. Again this feature does not depend 
significantly on the strength of the interaction parameters. Varying them in a broad 
range (factors between 1/5 and 2 for instance) leaves the histogram in Fig. 2.c es- 
sentially unchanged. This is another confirmation that the fluctuations, which we 
are looking at, reflect the generic properties of the system. 

The discussion presented above is based on the subspace of 2p2h states coupled 
to J n = + . Similar behavior can, however, be expected for the other multipoles. 
In smaller model spaces where the 2p2h diagonalization can be done reliably for 
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any multipole we see no qualitative difference when performing a similar study. Of 
course, in the smaller spaces, the mean square deviations are larger due to poorer 
statistics. 

Because of the similarity of wave functions, the most important doorway states 
for decay of a typical giant resonance are those 2p2h states which are generated 
by the single-particle basis of two major mean-field shells on both sides of the 
Fermi surface. Those states have, therefore, been discussed extensively above. To 
conclude about the degree of genuine chaoticity of the nuclear hamiltonian projected 
onto the full 2p2h space one needs, however, a more restrictive test. The point is 
that including further single-particle shells in the diagonalization may influence 
the higher-energy part of the spectrum studied above. Since such an extension 
cannot be made in practice we go the opposite way. Out of the total of 2696 
J 71 " = o + states we select a sequence of 400 states starting from 51-st state. Such a 
sequence covers the excitation energies up to about 30 MeV. Concerning the overlap 
with the high-energy 2p2h states, not taken into account in the diagonalization, we 
find this choice safe. The resulting NNS distribution and the A3 statistics [18] 
the latter being a measure of the rigidity of the spectrum are shown in Fig. 3, 
including all diagrams of Fig. 1. The A^(L) statistics is calculated as an average 
of the mean-square deviation from the straight line of the integrated density of 
(unfolded) states in the intervals of length L. A comparison with the corresponding 
Poisson and GOE predictions []|] is also made. We have verified that using an 
analogous sequence of states on the unperturbed level reveals a similar clustering as 
in Fig. 2. a. Furthermore, including diagrams (a) or (c) or both leads to decorrelation 
seen already for a string of eigenvalues of the same length. The presence of such 
a transition means that the final good agreement with the GOE does not involve 
any kinematical repulsion [19], but is entirely due to the dynamical correlations 
violating integrability and thus generating chaos. 

In conclusion, already on the 2p2h level the spectral fluctuations of the GOE 
can manifest themselves, provided the residual interaction is taken into account. 
Most important in this connection are the particle-hole rescattering effects which 
correlate the states from both sides of the Fermi surface. While the unperturbed 
2p2h states certainly provide a reasonable first approximation for the global energy 
location of the doorway states, they fail completely in the sense of fluctuations. 
This may have an influence on the fine structure of a resonance. The above aspect 
of nuclear dynamics should also be kept in mind when addressing the question (2(J 
as to whether the nuclear collective motion is Markovian or not pij| . In classical 
terms a chaotic system looses memory very fast because of exponential instabilities. 
It is thus natural to expect that chaos suppresses the role of the history also on 
the quantum level. Actually, the integration kernel in eq. ^ involves a sum of 
exponents eq. (0) over the whole background spectrum. For a spectrum with a 
strong tendency to clustering many terms in such a sum may add up constructively 
and thus amplify the r dependence in eq. @. The opposite should apply to the 
GOE spectrum. In fact, a simple estimate of such effects based on the spectra 
discussed in this letter shows one order of magnitude reduction of the amplitude of 
oscillations when going from the case of Fig. 2. a to Fig. 2.c. This is an important 
problem |22| which demands a more systematic study. 
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Figure Captions 



Fig. 1 : Diagrammatic representation of the two-body matrix elements in the space of 
2p2h excitations with explicit indication of the angular momentum coupling scheme. 

Fig. 2 : Nearest-neighbor spacing distributions (histograms) for the sequence of 
2p2h states coupled to J 71 = + as a function of the normalized relative distance 
s. The distribution has been generated in a basis of four major shells (two below 
the Fermi surface and two above) for a Woods-Saxon potential with parameters 
corresponding to 40 Ca. Part (a) displays the unperturbed case, (b) corresponds to 
the results from a diagonalization including diagrams (a) and (c) of Fig. 1 and (c) 
includes all the diagrams. The residual interaction used is given by eq. @. The 
dotted lines represent the Poisson and the solid lines the Wigner distributions. 

Fig. 3 : Nearest-neighbor spacing distribution (histogram in the upper part) and the 
A3 statistics (diamonds in the lower part) for a sequence of 400 low-energy states 
(between 51 and 450) obtained from a diagonalization of the residual interaction 
eq. @ in the basis of 2p2h states for J n = + . The dotted lines refer to the 
Poissonian spectrum and the solid lines to GOE predictions. 
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